The effects of fatigue on linear and angular kinematics during bilateral squat

This study was aimed to analyze in detail how the fatigue effects to kinematic parameters of body weight squat exercise (BSQ) by dividing a squat cycle into four different regions. Twenty-one male athletes participated in this study. Participants were divided into two groups according to their lower limb muscle ratio (LLMR). The BSQ was performed until participants were unable to continue the exercise due to the fatigue. Linear and angular kinematics were obtained by motion analysis software which has high validity and reliability. There was no significant but had large effect size interaction between fatigue conditions and LLMR groups in terms of knee ROM in the extension phase and hip angular velocity in braking phase of the flexion (0.08 > p >0.05, 0.18 > ηρ2 > 0.16). Fatigue condition did not have a significant effect on the duration in the acceleration and braking phases of BSQ (p > 0.05). There were many significant main effects on kinematics in the different regions due to the fatigue (0.01 < p <0.05, 0.44 > ηρ2 > 0.14). In the fatigue condition, there was a polynomial relationship between velocity of shoulder and hip joints (R2flex = 0.82, R2ext = 0.72) rather than linear (R2flex = 0.64, R2ext = 0.53) and coefficient correlations also decreased (rflex = 0.88 to 0.80, rext = 0.92 to 0.73). The sticking region was observed in the non-fatigue condition and disappeared when fatigue occurred. These results suggest that LLMR may be taken into consideration in the squat exercises, joint tracking may vary for velocity-based squat training and pre-test for sticking region observation may be apply with the BSQ.


Introduction
The bilateral squat is one of the most common exercises of developing strength [1,2].It also plays a vital role in the completion of everyday activities.Therefore, it is one of the most preferred exercises especially for the general population, beginners and athletes to rehabilitation and athletic development [3].It can be a form of rehabilitation as well [4].As squat is dynamic and closed chain, it is also influenced by a multitude of factors.The depth of the movement, positions of the feet on the ground, the strength of core muscles, the position of the head may alter the activation of muscles [2,[5][6][7].These are important factors that vary according to the individuals' goals.Additionally, men and women exhibit different motor strategies across the kinematic chain [8].
Fatigue is another critical factor in squat exercises [5,9,10].It has been shown in many studies that there are changes in the kinematics and kinetics of the movement due to the fatigue.Some of these; more forward leaning [6], forward shifting of the ground reaction force [11], reduction of range of motions [3] and decelerate the linear velocity of the joints [12].In resistance exercises, a phenomenon occurs in the fatigue condition.This phenomenon is called the "stick region".This region may be simplest summarized as the short-term slowdown and followed by acceleration during the lift (in the ascent phase).The method used in the studies is generally focused on the concentric phase related to this region.Because, the unsuccessful end of the movement or the observation of this "region" occurs at this stage [13].Even though much have been investigated so far, it is still not possible to say the exact reasons why this phenomenon occurred.In addition to the these, it is not possible to observe this region for each participant in the fatigue condition [13,14].Also, in this region, the acceleration and braking durations in the ascent phases are generally not equal and it is observed that the deceleration duration increases more with light external loads [15,16].Previous studies haven't emphasized the body weight squat exercise (BSQ) in exploring these phenomena.However, based on existing findings, it could be expected that the acceleration and braking durations would become more similar in this exercise.
Another important training effect related to squat is velocity-based training.Regarding the velocity of the joints, monitoring the bar speed and predicting one-repetition maximum (1-RM) methods have become very important and popular for the velocity-based training.However, BSQ was not investigated as much as weighted squats (i.e., back squat).Therefore, questions such as how the velocity change of the joints (i.e., hip, shoulder) during the exercise and whether the relationship between them is linearity or not, has not been fully answered.Another situation that needs to be considered here is that monitoring to the bar speed for upper extremity exercises can be more direct method (i.e., bench press).However, while the difference of joint velocities related to this situation has been proven [17,18], following the lower extremity development or training from monitoring the bar speed may be evaluated indirect method for a multi-joint movement (i.e., squat).In addition, lower limb muscle ratio (LLMR) has an important role in activities such as squatting, jumping, sprinting etc.Although the issue of how the LLMR develops with the different type squat exercises have been well researched, studies of how the LLMR affects the kinematics while the squatting are very limited.Given these gaps, investigating i) joint velocity kinematics during BSQ, ii) kinematic changes in participants with different LLMRs, and iii) comparing these kinematics in nonfatigue and fatigue conditions, can provide insights.The null hypothesis of this study, there was no tend to significant difference between the LLMR groups in terms of linear and angular kinematics.Therefore, the purpose of this study was to analyze in detail how the fatigue effects to kinematic parameters of BSQ by dividing a squat cycle into four different regions.

Participants
There were more men than women who met the inclusion criteria of this study.Therefore, all participants were selected as men to maintain group homogeneity.Twenty-one recreationally male athletes (age: 23.7±3.2 years (mean ± SD), height: 179.6±7.7 cm, body mass: 79.52±11.78kg) volunteered and gave their informed consent to take part in the study.All participants with a history of resistance exercise for at least two years (2 or 3 times per week).The inclusion criteria for this study were defined as (a) had already practiced and continued squat exercises in a resistance training program, (b) had no chronic limb discomfort or limitations that had a negative effect on BSQ kinematics (e.g., pes planus).The exclusion criteria defined as experiencing any injury during the BSQ and not being able to perform BSQ until exhausted.Participants in the study were subdivided into two groups according to the lean muscle mass of the lower limb.These data obtained from a body composition analysis device (Tanita BC 418, Tanita Corp., Tokyo, Japan) which is using the bioelectric impedance method in order to reveal the effect of lower limb muscle mass ratio in the fatigue condition.Based on this, the lean muscle mass (lower limb) data were divided into two groups as low and high by making cut off at 50%.The range of the low group (n low = 10) was 32.00%-33.86%and mean = 33.28%and the range of high group (n high = 11) was 34.25%-36.63%and mean = 34.85%.The range of the whole group (n = 21) was 4.62%.The experimental procedures undertaken were approved by Ege University's Medical Research Ethics Committee and are in agreement with the principles of the Declaration of Helsinki.The participants were briefed on the study procedures and informed consent was obtained before they took part in the study.

Instrumentation and data processing
The BSQ was performed using a tempo standardized to a two-count cadence at 45 beats per minute (bpm).The cadence was 45 bpm (1:0:1:0) which signifies a one-count eccentric motion followed immediately by a one-count concentric motion, without any count for rest before beginning the next repetition.Linear and angular kinematics of the shoulder, hip and knee joints were examined during the performance.Escamilla et al., (2001) previously described as i) descent phase, ii) acceleration phase, iii) sticking region, iv) maximum strength region and v) deceleration phase for a squat cycle with external load.However, since this study is performed with only body weight, a squat cycle divided into four different regions.These regions were generated in order to make a body weight squat-specific assessment (described below, Fig 1).Another important reason is that the sticking region is not expected to be observed in BSQ as mentioned above.To assessment for kinematic effects of BSQ; the third repetition after the initiation and the previous one before the last accomplished repetition were analyzed in order to observe the movement pattern properly (total; two).The ROMs (˚), angular velocities (rad.s - ), vertical velocities (m.s -1 ), vertical accelerations (m.s -2 ) and time (s) analyzed for each region.Vertical linear and angular velocities (m.s -1 , rad.s -1 ), ROMs (˚) and accelerations (m.s -2 ) except peak velocities, angles and time were averaged.The first position of upward position (neutral standing posture, start-position) took as initial angle via motion analysis software (full extension, initial angle for hip and knee:180˚, ankle:90˚), the last point of downward position (descent) evaluated as flexion phase and from the last point of downward position to the last point of upward position evaluated as extension phase.The differences between these angles allowed us to obtain the ROM of hip, knee [180˚minus the last point of downward angle (V flex region, ROM flex ); the last point of downward position angle minus 180˚(V ext region, ROM ext ); angles between; acromion-trochanter major for hip, trochanter major-malleolus lateralis for knee, epicondylus lateralis-(malleolus lateralis)-malleolus lateralis for ankle].This study was carried out in a laboratory in two sessions at 72-hour intervals.In the first session, after the participants had done a 10-minute warm up, retroreflective markers were positioned for the acromion (shoulder), trochanter major (hip), epicondylus lateralis (knee), malleolus lateralis (ankle) and caput metatarsalis (toe) points before the exercise.The BSQ was performed for adapt to the metronome and markers.Those who could not adapt to the metronome were excluded from the study.In the second laboratory session, retroreflective markers were placed at the same reference points as in the first laboratory session.The BSQ was captured by a Basler Ace acA1300-200uc camera (Basler AG, Ahrensburg, Germany) with a frequency of 100 Hz (0.01 s).A digital metronome was used to standardize the participants' BSQ pace (45 bpm, 1:0:1:0).The camera was placed on the right lateral side of the participants (dominant leg), distance: ~2m, camera angle according to the participant: 90˚ [19].After the participants had enough warmed up, they performed the BSQ looking straight ahead at a fixed point that had been placed on the wall with their feet approximately opened hip width (greater trochanter width) with feet externally rotated to a maximum of 15˚.Because, when comparing 2-D and 3-D analysis for the squat exercises, foot angle (i.e., 0-15˚) and narrow stance (107±10% of shoulder width) form of squat are showed minimum error between 2-D and 3-D analysis (Escamilla et al., 2001).For this reason, in this study, squat performed that can be called "narrow stance" was carried out.The participants raised their arms at the descent phase and released their arms the arms during the ascent phase to avoid the fatigue of upper limb.The flexion phase of the squat was not limited.We identified fatigue as a situation in which participan's joint kinematics are severely disrupted, failed to keep up with the metronome and cannot continue the exercise (unable to squat or stand up).

Defining the regions for BSQ
In order to examine each phase of the BSQ more details, a repetition is divided into four regions (Fig 1).The first region was from the beginning of flexion (full extension) to the peak velocity of flexion (V 1down ).Acceleration values of first region was also called "Descent Acceleration".The second region (V 2down ) was from peak velocity of flexion (V downpeak ) to the end of flexion.Acceleration of second region corresponding to the same phase denotes "Descent Deceleration ".The third region (V 1up ) started from beginning of extension to peak extension velocity (V uppeak ), acceleration values of third region was also called "Ascent Acceleration".The last region was from peak velocity of extension (V uppeak ) to end of extension which was also called V 2up (Fig 1).The averaged acceleration value corresponding to the same phase denotes "Ascent Deceleration".All velocity (m.s -1 , rad.s -1 ) and accelerations (m.s -2 ) except peak velocities, angles and time were averaged.

Statistical analysis
The present investigation was designed as a cross-sectional, between-within subject study.The sample size was calculated using GPower 3.1.9.7 [power size (1-β) = 0.93, the effect size (f) = 0.40, type-1 error (α) = 0.05, number of groups = 2, number of measurements = 2, cor.among rep.measure = 0.5, correction ε = 1, input/direct parameter Z 2 r = 0.14].All negative variables are shown as positive by taking absolute values.To assess differences between the 2 squat conditions (non-fatigue, fatigue) x 2 LLMR group (low-high), analyses of variance (twoway mixed ANOVA) was performed.Greenhouse-Geisser adjustments of the p-values were reported due to sphericity assumption was violated [24].Partial eta squared (Z 2 r ) was used as an effect size and categorised as small (0.01-0.06), medium (0.06-0.14) and large effect (> 0.14) [25].Test-retest reliability for the linear and angular kinematics between the nonfatigue and fatigue conditions was analysed by ICC.The ICC are reported with 95% of confidence interval (95%CI) and categorized as poor (< 0.5), moderate (0.5-0.75), good (0.75-0.9) and excellent (> 0.9) [26].The linear and polynomial regressions were performed to establish the velocity relationship between shoulder and hip joints in fatigue conditions.Pearson correlation coefficients (r) were also used to determine the significance of the velocity of shoulder and hip joints relationships.Statistical analysis and graphics were performed by using RStudio 1.3 [27] with five packages [28][29][30][31][32].The regression graphics created with Microsoft Excel (v16.0,Microsoft Corporation, WA, USA).The statistical significance level was set at p � .05.

Results
In this study's analysis; data were screened for assumptions and outliers.Out of twenty-six significant data, nine were outlier as assessed by boxplot.When the results were compared with the previous results (with-outlier), only small differences were observed.Nevertheless, corrected data (without-outlier) were reported.There were also eight non-significant (.08 > p > .05)but large effect size variable (Z 2 r > .14).These variables are also included in the report (Figs 2C, 3D-3E, 4C and 5D-5F).There was homogeneity of variances (p > .05)and covariances (p > .05),as assessed by Levene's test of homogeneity of variances and Box's M test, respectively.Participants performed 72 ± 27 repetitions during the whole exercise.However, only two repetitions were taken into account as described above (non-fatigue vs fatigue).
The main effect of all measurements' mean and standard deviation are shown in Table 2.The statistically significant main effects (n = 21, non-fatigue vs fatigue condition) of linear and angular kinematic ANOVA results (above the box-whisker plots) are also shown in Figs 2-5.
The test-retest reliability for linear kinematics of shoulder, hip and knee joints between the non-fatigue and fatigue conditions found that there was poor to good reliability (0.26 < ICC < 0.80) and there was poor to good reliability (0.08 < ICC < 0.79) for the angular kinematics of hip and knee joints between the fatigue conditions (Table 3).
The main effect of exercise duration did not show statistically significant difference in the four regions (and in the flexion and extension phases) of the exercise at the different fatigue condition (p > .05)(Table 4). .The linear velocity of the shoulder and hip joints was drawn and examined for each participant.Meanwhile, it was found that one of the participants draw a curve similar to "sticking region", and even in the non-fatigue condition (Fig 8), when the video analysis of this participant is examined carefully; it took exactly 0.06 s in the pre-sticking point (from V 0 to V max1 ).Then, the time between V max1 and V min were observed as 0.1 s in the sticking region.The ROM of the participant at the end of the flexion phase (V 0 ); hip: 72˚, knee: 67˚, ankle: 67˚, and the participant reached from V max1 to V min ; hip: 72˚, knee: 70˚, ankle: 70˚(Fig 8).There were also found similar curves for the shoulder joint as in shown Fig 9.

Discussion
In this study, we investigated linear and angular kinematics of body weight squat exercise in two different groups according to their lower limb mass ratio (low and high) between the nonfatigue and fatigue conditions.They performed the exercise until exhaustion occurred.We showed that different lower limb mass ratio effects on kinematics of exercise.There was not observed any statistically significant duration difference at the different regions of exercise while there was significant different in the kinematics.In the fatigue condition, we showed that there was polynomial relationship between shoulder and hip joints.We posit that this could provide a new approach to velocity-based training.Additionally, our observations suggested that an individual's habitual movement form during the squat exercise might contribute to the occurrence of sticking regions.
The mean flexion velocity was greater in the low group than the high group.It may be considered that the low group is trying to compensate the fatigue to comply with the tempo by increasing ROM during the extension phase.On the other hand, although the same trend was observed for the hip, no significant interaction was found, because the increases were less (ROM_Ext).Both groups are also tried to compensate for increased the initial phase of angular velocity (V 1down _Deg_Velo).This suggests that individuals with high lower limb mass ratio (LLMR) are more likely to be able to control more in the flexion phase during the body weight squat exercise (BSQ).It may be concluded that more caution should be taken with individual's LLMR below 30% for BSQ.

The main effects of angular and linear kinematics (n = 21, non-fatigue vs. fatigue condition)
Fatigue resulted in a decrease in the angle of peak velocity (V downpeak _Angle) for both the hip and knee joints during the flexion phase.This indicates that in the fatigue condition, the hip reached peak velocity closer to the ground in this phase.This decreasing in the hip and knee is due to the fact that the participants performed the flexion phase (V 1down _Deg_Velo) with faster angular velocity (ω).In the flexion phase of fatigue condition, there was increasing in the angular velocities of the hip and knee joints in the initial acceleration region (V 1down _Deg_-Velo) as well.However, during the deceleration phase (V 2down _Deg_Velo), the angular velocity of the hip decreased while the knee joint's angular velocity increased.This may be due to the different mobility characteristic of the hip and knee joints.
In the non-fatigue condition, participants were able to perform a rapid flexion followed by a slower extension phase to keep up with the exercise pace.However, it resulted in an increase In the "A-B" represent peak angles were observed.In the "C", it represents the ROM of knee joint.In the "D-E-F", they represent mean velocity of knee joint in the regions.There was statistically significant main effect of the velocities in the "B-C" (p � .05,Z 2 r > .14).There was not significant but had large effect size in the "A-D-E-F" (p > .05,Z 2 r > .14).https://doi.org/10.1371/journal.pone.0289089.g005 of 4.62% (flexion) and 23.58% (extension) in the fatigue condition.On the other hand, it is observed that mean flexion velocities (V flex ) were less than the extension velocities (V ext ).This condition particularly supports the effectiveness of fatigue during the flexion phase.Although the contraction phase (extension) is a general concern, especially in the studies related to velocity, this situation which increases even with the body weight, supports that should be more cautious in the flexion phase (descent) with heavy loads.In these conditions, it is usual to observe the increasing the load on the joints [3,5,33].This situation may cause undesirable results in long term (i.e., injury).In the fatigue condition of the extension phase, decreased differences in the angles of peak velocity (V uppeak _Angle) was observed in the hip and knee joints.This occurred when the angular velocity (V 1up _Deg_Velo) of the hip and knee joints increased in the extension phase.This increased velocity caused the peak velocity to be reached earlier in the extension phase.However, when the exercise durations were examined, neither acceleration (V 1down(flex.)-V 1up(ext.) ; NF vs. F) nor braking (V 2down(flex.)-V2up(ext.) ; NF vs. F) phase did not show difference.In some studies, in order to avoid the effect of the stretch and shortening cycle at the end of the flexion phase (between eccentric and concentric phases), a short pause period is implemented approximately between 2-4 s [15,34].No such method was followed in the present study.Therefore, the observed increases in angular velocity may be attributed to the stretch-shortening cycle [35,36] or the participants' efforts to synchronize with the exercise pace.Another factor to consider is inertia.The contribution of the velocity in the extension phase to the increased in moment may be explained by faster flexion phase occurred (e.g., rebound effect of the movement).Hence, introducing a pausing period at the conclusion of the flexion phase could potentially yield positive effects, akin to isoinertial strength training, known for its performance-enhancing benefits [10,37,38].The present study has been one of the few studies that observed to higher vertical linear and angular velocities in the fatigue condition.There were also increased difference in the flexion ROM of the hip and knee joints in the fatigue condition.This may be considered that participant try to keep themselves upright during the descent phase and release themselves faster due to the effects of gravity and fatigue.Another situation to consider was that in case of fatigue, the ROM of the knee joint would decrease during the descending phase, and the trunk could be excessive forward lean [39].However, this possibility disappeared with the statistically significant increase in the flexion ROM of the knee joint without this increase in the hip joint.
In the flexion phase of the linear velocity, significant difference was observed in the shoulder, hip and knee joints.While this difference appears in linear velocities, no such difference was found in angular velocities (Flex_Deg_Velo) between fatigue conditions.However, flexion phase of the ROM supports this situation.While the linear velocity of the hip and knee showed difference in the extension phase, this was not observed in the shoulder joint.This suggests that participants exerted effort to maintain an upright posture during the fatigue condition.
Otherwise, significant difference in this linear velocity could be observed during further forward and backward movements (i.e., excessive forward lean).Taken together, the non-significant difference in the ROM of hip supports this situation.The velocity relationship between shoulder and hip joints in the fatigue condition In the studies of squat exercise related to velocity, it is generally performed by following the bar velocity [12,14,[40][41][42].The excessive forward lean of the trunk could be observed in the free weight back squat exercise [43], as a result of this, it might have a negative effect on the  condition.At the same time, the correlation coefficient between the velocities of the two joints in both conditions also decreased.In other words, fatigue condition in the extension phase caused further deterioration of the velocity relationship between the two joints.In fact, a study showed that monitoring bar speed for 1-RM with light loads did not show very strong linear equations as well [44], and Izquierdo et al. [45] found that parallel squat exercise had a significantly higher mean velocity in the last repetitions (60-75% 1-RM) than bench press exercise.
In the present study performed with only body weight, hip velocities were lower than shoulder velocities in the each of the flexion and extension phases.According to present study, there may be more logical way to monitor hip joint instead of shoulder (bar) for velocity-based trainings and especially predict 1-RM with low-loads when fatigue occurred.Therefore, velocity-ranges [46] might be revised again for squat exercise or the desired exercise intensity would be observed more accurately.The present study raised doubts about observing the upper extremity for the squatting exercise may be considered as a partially indirect method.Therefore, the differentiation of the monitored region (bar vs hip) may allow better comparison of these velocities.This may give a new perspective for further studies.The lack of strong R 2 (i.e., R 2 < .95)prevented that generating a generally acceptable and more predictable power of the equations in this study.However, in the cases where the number of sample size is higher, it seems likely to find stronger predictable equations with higher R 2 [40].Although it is more logical and reliable to monitor bar speed for 1-RM estimation methods with heavy loads [44,[47][48][49], this may not benefit the same rate in velocity-based squat trainings in general.The differences in these velocity outputs (hip and shoulder) can be considered quite normal that the joints take different distance at the same period of time.However, the situation that needs to be emphasized here was the R 2 and coefficient correlation (r) values which are changed in the different fatigue conditions.This situation, together with the external load, may cause less difference in the shoulder and hip joints.However, in an exercise with free weight, excessive movement of the trunk in the fatigue condition may prevent us from observing the desired velocity from shoulder (bar).According to this study, it is more appropriate and predictable for the relationship between shoulder and hip joints had non-linear relationship with a polynomial curve after fatigue occurred.

Is it possible to observe "sticking region" for BSQ?
The "sticking region" could be defined as deceleration and re-acceleration when lifting the bar with additional load in the contraction phase (moving upwards) during exercise (i.e., back squat, bench press) and usually observe before completing the movement, which may be result of exhaustion/inability to perform the exercise [13,40,50,51].In this study, it is observed that one of the participants draw a curve similar to "sticking region" with only body weight and this happened within just 3˚and 0.16 s (from V 0 to V min ).Therefore, it was not noticed by the researchers during the test that the participant had such a tendency.It is surprising that was observed without fatigue in the BSQ, although sticking region in the literature is generally occurred at the end of exercise (last repetition) with additional loads, and in some cases not even observed [13,14].Another interesting aspect is that during the flexion phase in the fatigue condition, a similar slowdown is observed while descending.This may be a situation-specific to the hip joint and this is observed in the graphics of some studies.Kristiansen et al. [52] showed that there is no generally acceptable optimal lifting method, each individual may find their own way to successful lifting and individual differences may affect performance positively such as shorter trunk and thigh length.In this study, the participant whose curve is drawn, did not have any disability for proper form of the BSQ.It has been observed that both joint (shoulder and hip) draw similar curves.It may confirm the speculation about the monitoring of velocity in this study.Although they had similar tendencies, as a result of a detailed examination, a linear relationship was not observed between two joints.Also, observing the sticking region in the nonfatigue condition of BSQ, it may arise the possibility that this phenomenon may occur participant's form of squat (unconsciously happening, performing exercise in their own way).This situation, which corresponds to 4.76% of the participants in the present study, considering that it is performed with external load, may increase the possibility of observing the sticking region.R. van den Tillaar et al. [53] found that this region observed in two-thirds of the participants in 6-RM squat exercise.Martı ´nez-Cava et al. [54] found that sticking region was not occurred in the half squat even with heavy loads.Perhaps before testing this phenomenon with external loads, the tendency of participants to such situation may be compared by performing a pre-test examination with the BSQ.

Practical applications
The lower limb muscle ratio differences may lead to significant interactions for angular kinematics in the fatigue condition.If individuals have less than 30% LLMR may increase the risk of movement controlling especially in the flexion phase after fatigue occurred.In this sense, it may be important to calculate the LLMR before squat exercises.The relationship between the shoulder and hip joints may be better explained by polynomial rather than linear relationship in the fatigue condition.Because of this situation, tracking the hip joint may provide an important recommendation for observing squat exercise in the velocity-based training.The studies with higher number of participants may contribute more for the future researches about this topic.Depends on the present study, it may be necessary to pre-test with BSQ before studies for observing the sticking region.Therefore, individuals who tend to show sticking region may be identified and the studies carry out with this method may provide more contribution as to why this phenomenon occurs.

Fig 2 .
Fig 2. Main effect of the linear velocities of the shoulder joint.ANOVA results showed upon each figure.Black lines of box plots in these figures represent medians.The numbers inside the boxplots represent the mean velocities.There was significant main effect of velocities in the "A-B-D-E-F" (p � .05,Z 2 r > .14).There was not significant but had large effect size in the "C" (p > .05,Z 2 r > .14).https://doi.org/10.1371/journal.pone.0289089.g002

Fig 3 .
Fig 3. Main effect of the linear velocities of the hip joint.ANOVA results showed upon each figure.Black lines of box plots in these figures represent medians.The numbers inside the boxplots represent the mean velocities.There was statistically significant main effect of the velocities in the "A-B-C-F-G" (p � .05,Z 2 r > .14).There was not significant but had large effect size in the "D-E" (p > .05,Z 2 r > .14).https://doi.org/10.1371/journal.pone.0289089.g003

Fig 4 .
Fig 4. Main effect of the linear velocities of the knee joint.ANOVA results showed upon each figure.Black lines of box plots in these figures represent medians.The numbers inside the boxplots represent the mean velocities.There was statistically significant main effect of the velocities in the "A-B" (p � .05,Z 2 r

Fig 5 .
Fig 5. Main effect of the angular velocities of the hip and knee joints.ANOVA results showed upon each figure.Black lines of box plots in these figures represent medians.In the "A-B" represent peak angles were observed.In the "C", it represents the ROM of knee joint.In the "D-E-F", they represent mean velocity of knee joint in the regions.There was statistically significant main effect of the velocities in the "B-C" (p � .05,Z 2 r > .14).There was not significant but had large effect size in the "A-D-E-F" (p > .05,Z 2 r > .14).https://doi.org/10.1371/journal.pone.0289089.g005

Fig 6 .
Fig 6.The relationship between the hip and shoulder joints in the different fatigue condition of the flexion phase.https://doi.org/10.1371/journal.pone.0289089.g006

Fig 7 .
Fig 7. The relationship between the hip and shoulder joints in the different fatigue condition of the extension phase.https://doi.org/10.1371/journal.pone.0289089.g007